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Abstract 

We study chiral symmetry breaking in QED3 with Nf flavors of four- 
component fermions. A closed system of Schwinger-Dyson equations for 
fermion and photon propagators and the full fermion-photon vertex is pro- 
posed, which is consistent with the Ward-Takahashi identity. A simplified 
version of that set of equations is reduced (in nonlocal gauge) to the equation 
for a dynamical fermion mass function, where the one-loop vacuum polariza- 
tion with dynamically massive fermions has been taken into account. The 
linearized equation for the fermion mass function is analyzed in real space. 
The analytical solution is compared with the results of numerical calculations 
of the nonlinear integral equation in momentum space. 

ll.30.Qc, ll.30.Rd, ll.lO.Kk, ll.15.Tk 
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Typeset using REVT^ 



I. INTRODUCTION 



Quantum electrodynamics in one temporal, and two spatial dimensions with Nf fla- 
vors of four-component Dirac fermions (QED3) continues to attract attention as a useful 
field-theoretical model for studying such phenomena as confinement and chiral symmetry 
breaking (xSB) which are out of the realm of perturbation theory. The model has properties 
reminiscent of QCD and other four-dimensional gauge theories. Thus it has an intrinsic di- 
mensionful parameter e^, the coupling constant, that plays a role similar to the QCD scale A, 
and the efTective coupling a{q) approaches zero at large momenta q [1]. Studying Schwinger- 
Dyson equation (SDE) for the fermion self-energy under the bare vertex approximation in 
the Landau gauge, Appelquist, Nash and Wijewardhana [2] have shown the existence of 
a finite critical number of fermions, Nc, below which the chiral symmetry is broken and 
fermions acquire a dynamical mass. In the vicinity of Nc this mass is exponentially small 
with respect to the natural dimensionful parameter e^, showing how a hierarchy of scales 
can occur in gauge theories. These results were confirmed by Dagotto, Kogut, Kocic [3,4] 
in their Monte-Carlo simulations of non-compact lattice QED3. 

Continuum studies of x^B use SDE which, being an infinite chain of equations, must 
be truncated in some way. The bare vertex approximation was criticised [5,6] as not being 
consistent with the Ward-Takahashi (WT) identity. An effective way to truncate SDE is to 
make an ansatz for the vertex satisfying certain criteria [7] (for a review of SDE in QED3 
see [8]). While at the present time the latter approach seems to be the most effective way to 
satisfy WT identity under truncation, it does not shed much light on an approximation from 
the physical point of view. Ideally, one should solve SDE for the vertex itself, but that is 
quite a formidable task. However, in the second section we propose one possible truncation 
scheme consistent with WT identity not mentioned before in the literature. 

Another important point in studying the Nf dependence of x^B is that it is necessary 
to insert vacuum polarization effects into the equation for the fermion dynamical mass (in 
quenched approximation, chiral symmetry is broken for any Nf [9,10,4]). In expansion 
the one-loop vacuum polarization with massless fermions is used. However, as was men- 
tioned in [11,12], the inclusion of massive fermions changes drastically the infrared behavior 
of the model, leading to the appearance of a confining (logarithmic) term in the potential at 
large distances. To take into account this confining nature of the model one can study the 
coupled system of SDE for the fermion self-energy, the photon polarization tensor and the 
vertex (for example, to study the coupled system of equations for them mentioned above). 
In this paper we chose another approach, replacing the fermion dynamical mass function 
in the polarization tensor by a constant dynamical mass of a fermion. While keeping the 
confinement property of the model this makes analysis much simpler. Because this approx- 
imation just corresponds to using the bare vertex, in order to overcome the inconsistency 
with the WT identity (or, at least, minimize its violation) we selected a nonlocal gauge 
[13,14] to keep the fermion wave function Z{p) = 1. 

In the third section we reconsider the analysis of linearized SDE for a fermion mass 
function, converting it into a Schrodinger-like equation in real Euclidean space with an 
effective potential. For the vacuum polarization with massless fermions this potential be- 
haves as — A/r^ at large r (A = 32/37T^Nf). From the theory of the Schrodinger equation 
for such kinds of potentials, it is known that it has an infinite number of bound states if 
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A > l/A[Nf < Nc = 128/37r^). However a finite number of bound states might exist for 
A < 1/4 as well. We prove that for the potential under consideration there are no bound 
states for A < 1/4, which means that the critical number Nc does indeed correspond to the 
appearance of an infinite number of bound states. 

In the case of massive fermions, the vacuum polarization leads to a Coulomb tail in 
an effective potential at large r, however, as we show in Section 4, this does not alter the 
exponential-type behavior of a dynamical mass near the critical Nc- This is, of course, 
in agreement with bifurcation analysis (see, for example, [15,16]). The main difference 
between a vacuum polarization with massive or with massless fermions lies in an overall 
scaling factor, which is larger for massive fermion loops in vacuum polarization. Prom the 
physical point of view it means that chiral symmetry breaking takes place essentially at 
intermediate distances 1/a < r < 1/m {a = e'^Nf/8) and is not influenced much by large 
distances, where the confinement property of the model exhibits itself. This situation is 
similar to that in QCD4 [17,18], where the chiral symmetry breaking scale is also different 
from the scale of confinement. 

The real space approach allows us to get a good analytical solution, fitting nicely com- 
puter calculations when Nf is close to Nc. In Section 5 we return to momentum space to 
analyze nonlinear SDE for a fermion dynamical mass function. We give the results of com- 
puter calculations of the dynamical mass function, which show very convincingly how a new 
energy scale (mass of a fermion) appears that is much smaller than the natural scale (e^) of 
the model. This might be important for a better understanding of a hierarchy of scales in 
unified theories. 



II. TRUNCATION OF SCH WINGER-DYSON EQUATIONS 



In QED3 the Schwinger-Dyson equation for the fermion propagator in Minkowski space 
is given by 

^-'(p) = P- 7^ / d'krS{k+p)r^(k+p,p)D^,(k), (1) 
and the SDE for the photon propagator is 

^;J(9) = ^vJ(g) + n^.(9), (2) 

where 

• 2 

^'""{q) = 7^ / Tr [7^^(fc)r (fc, k - q)S{k - q)] . (3) 



M 



In Eqs.(l) and (3) T'^{k,p) is the dressed fermion-photon vertex which should satisfy the 
Ward-Takahashi identity 

{k-p)^r^{k,p) = s-\k)-s-\p) (4) 

Eqs.(l) and (3) for the full fermion and photon propagators should be solved together with 
the equation for the full vertex function r^{k,p) 
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which in its turn contains an unknown kernel K{p, k, q) (fermion-antifermion scattering 
amphtude), and so on. Because the system of SDE is in fact an infinite chain of integral 
equations for n-point Green's functions, we are forced to use some truncation scheme, which 

we require to be consistent with WT identity. One can verify that if we take the bare 
vertex approximation (F^ = 7^^) in Eq.(l) and use the ladder approximation for the kernel 
K{p,k,q): 

K{p + q,k + q, q)cd,ba = 7ad7c6^Ap(g), (6) 

we get a closed system of SDE consistent with WT identity (4) (Eqs.(2),(3) remain unal- 
tered), see Fig. 1. Indeed, let us multiply Eq.(5) with the kernel (6) by {k —p)n, then taking 
into account Eq.(l) (with — 7^) we write 



k-p = S-\k) S-\p) 
J d'qrS{q + k)YD,.{q) -7^/ d'q^Siq + p)Y (7) 



^ ^ M ^ ^ M 



{k - p),r^{k,p) - S-\k) + S-\p) = 7^ / d'qi'S{k + q) 



which allows us to reduce Eq.(5) to the following: 

X \k - p)^j:^{k + q,p + q)- S-\k + q) + S-\p + q) 

xS{p + q)YDM- (8) 

It is evident now that Eq.(8) has a solution satisfying WT identity Eq.(4). 

Due to the Ward-Takahashi identity, the vacuum polarization tensor Eq.(3) has the form 

= i-9f.uq^ + q^QuMq^), (9) 

where n(g^) is the vacuum polarization (we recall that we use four-component spinors). 
The general form for the full fermion propagator can be expressed as 

where Z is the fermion wave function and M is the mass function. Then the equation for 
the mass function is 



^(P') ie"" „ ,M(k^)Z(k^) , , 



M 

and for the fermion wave function 
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p 1- 



1 



) 



(27r) 



,2 



M 



(12) 



The general tensor structure of the vertex contains eight scalar functions [7] ; thus, together 
with the functions Z{p^), M{p^) and n(p^), we can write down a coupled system of integral 
equations for eleven unknown scalar functions. Perhaps, this is the simplest truncated set 
of SDE consistent with the WT identity. We emphasize that such a truncation scheme can 
be used also in four-dimensional QED with one (though important) difference: due to the 
superficial linear divergence of the electron self-energy we should take a photon momentum 
as an integration variable in order to avoid shifting variables. Certainly, this scheme does 
not satisfy the requirement of multiplicative renormalizability, but that property is a feature 
of the whole theory rather than of any approximation scheme. 

Now, it is well known that Z is a gauge-dependent function, and we can use the freedom 
of choosing a convenient gauge to make the function Z{p'^) — 1. To achieve this we use 
a nonlocal gauge fixing procedure [13,14]. The photon propagator in a general nonlocal 
covariant gauge is defined as [19] 



with ^{q'^) being a function of momentum q rather than a constant. Simmons proved [13] 
that a suitable form for ^(g^) is the following: 



With this gauge function C{q'^), the right-hand side of Eq.(12) vanishes when averaged over 
the direction of k. 

For Eq.(ll) we obtain, in Euclidean formulation. 



where q^^ = k^^ — p^. 

The system of Eqs.(2), (3), (5), (6) and (15) is still far too complicated for an analytical 
study, and we postpone the full investigation of it to the future. Here we note only that the 
mass function M(p^) is connected to scalar functions from the vertex through the vacuum 
polarization n(g^). If we assume that this connection is not crucial, we can take the one- 
loop approximation for the vacuum polarization, replacing the vertex by the bare vertex 
7^^, and the running mass function MijP') by its value M(0). This approximation, though 
keeping the transversality condition for Il^y[q), leads to violation of the WT identity Eq.(4), 
however, in the nonlocal gauge the WT identity is still approximately satisfied. 

Thus the vacuum polarization in the one-loop approximation takes the form (with Eu- 
clidean momentum q): 




(13) 




(14) 




(15) 




(16) 
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where in our case we take m = M(p^ = 0). Taking massless fermion loops in vacuum 
polarization, we find 

n(g2)= « = (17) 

q 8 



III. COORDINATE SPACE FORMULATION 



First we study Eq.(15) in the linearized form, i.e. 



E 



^2 _|_ ^2 q2 



i + n(g2) 



(18) 



where we have replaced M{k'^) in the denominator by m = M(0). In terms of the coordinate 
space function 



d-V M(p2^ 



(19) 



the equivalent of Eq.(18) is 



d^pd^^A; M(A;2) 1 



1 + n(g2) 



+e(a 



(20) 



Performing a shift of integration variable, p ^ p + k, we find a Schrodinger-like equation for 
the function ■0(r): 



where the potential V{r) is defined as 



(21) 



y(r) = 



27r2r 



dp 



smpr 



1 + n(p2 



(22) 



In order to simplify (22), for both the massless and massive case we approximate the vacuum 
polarization by 



a 



q + P 



(23) 



where /S = am, a being some fitting constant, and a = Nfe'^/8. The massless vacuum 
polarization corresponds to putting a = 0, while the value a = Src/A corresponds to approx- 
imating both infrared and ultraviolet behaviors of the exact one-loop expression for 11 (g^). 
The above approximation allows us to perform the integral in ^{q^) analytically, i.e. with 
(23) 
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(24) 



With a — 1, i.e. = the potential can be expressed as 



/ 



dx 



sm 



X 



a; 



X + 1 + P X 

The last expression can be rewritten in a more convenient form: 
4/3 16 



1 + /?. 



(25) 



\/(r) = 



X 



1 



1 

+ ^ - 



dx sin 



1 



{1 + P)ix + 1 + /3) x^ 2{1 + P)x 



x-^ 



log 1 + 



1 + 3 



(26) 



Wc can get rid of the oscillatory behavior of the sine by performing a contour rotation in 
the lower right quadrant of the complex-plane of x. So 



V{r)-- 



/9 



16 



X 



1 



dy exp{—yr) 







1 + /3 



(27) 



Let us consider first the case of vacuum polarization by massless fermions, i.e. m = 0, which 
coincides precisely with the approximation used by Appelquist, Nash and Wijewardhana [2] 
in momentum space. The potential (27) takes the form 



16 7 

This potential behaves at large distances as 



y2 y2 



32 



1 

— arctan y 

y6 



(28) 



V(r) ~ -— Ai = 
whereas at small distances it has Coulomb-like behavior: 

r TcNf 



(29) 



(30) 



Prom the theory of the Schrodinger equation for such kinds of potential, it is known [20,21] 
that there is an infinite number of bound states if Ai > 1/4 [Nf < Nc = 128/37r^). However 
a finite number of bound states might exist for Ai < 1/4 as well. Let us show that for the 
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potential under consideration there are no bound states at all for Ai < 1/4. Indeed, one can 
prove the following inequality 



V{r) > / dy — > — J. 31) 

r J y + 1 







Since the Schrodinger equation 

l^-V^ - ip{r) = -mV(r) (32) 

has no bound states for Xi < ^ [20], it is evident from the inequality (31) that the Schrodinger 
equation (21) with the potential V{r) (28) has none either. Thus we have proved that the 
critical number of fermions, 

1 28 

A^c = 3^ ~ 4.32304, (33) 

does indeed correspond to the appearance of an infinite number of bound states in an 
effective Schrodinger equation in r-space, if we consider the vacuum polarization by massless 
fermions. The critical number of fermions (33) coincides with the result obtained by Nash 
in a different way [22] and that was claimed as being gauge invariant to leading order in the 
1/Nf expansion. 

The quasi-classical approach gives the following spectrum of bound states for potentials 
behaving at large distances as — A/r^ [21]: 



En = -ml ~ - exp 



27rn 



V^i - 1/4 



(34) 



at n >> 1, coinciding in fact with the result of Refs. [2,22]. As was pointed out in [2], this 
spectrum is closely connected with approximate scale invariance of QED3 at intermediate 
distances 1 < r < 1/m ^. However the quasi- classical theory does not tell us anything 
about lower levels of the Schrodinger equation. It is clear, also, that the proof given above 
breaks down in the case of massive fermion vacuum polarization, since the potential (27) 
has a Coulomb-like tail at large distances. The latter is simply a reflection of the fact that 
including a dressed fermion propagator restores the confinement property of QED3 in the 
sense that a logarithmic term reappears in the real potential at large r [11]. In the next 
section we shall solve the Schrodinger equation with the potential (27) and study the scaling 
properties of the dynamical mass near Nc- 



^We note that in QED4 it is exactly the scale invariance of the ladder and/or quenched approxima- 
tion which is responsible for exponential-like behavior of a dynamical fermion mass. For physical 
explanation of this fact see Ref. [23]. Further discussion of the role of scale invariance in QED4 
can be found in Refs. [24-26]. 
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IV. SCALING PROPERTIES NEAR CRITICALITY 



In this section we study the scahng properties of the dynamical mass near Nc, for both 
massless and massive vacuum polarization. For large r, V{r) (27) has the following form 

V{r)^-^-^, r»0{l), (35) 

where 

A3 = — — , p = am. (36) 

In deriving this we neglected terms /3 with respect to 1, since we consider V{r) near Nc, 
where m is very small (and so is (3). For large r we have thus a — potential with a small 
Coulomb potential which starts to contribute for r > For small r, V{r) behaves like a 
Coulomb potential plus a constant part 

V{r)^-^ + hi, r<C»(l). (37) 
r 2 

We split our problem into two regions, for large r, and small r. Neglecting angular depen- 
dence, which is not relevant for our purpose, we obtain 

-r^ + -—+— + m V'^O, r>l, 38) 

dr^ r dr y r y 

d^^ ^di) (\2 , ^ 

dr^ r dr y r y 

where = |Ai + w?. With ■^(r) = f{r)/r, the equations can be expressed as Whittaker 
equations 



dx^ \ 2a; 4^ 













2a; 










dx^ 


2fhx 



/ = 0, X = 2mr, r > 1, (40) 
/ = 0, X = 2mr, r < 1, (41) 



where the general form of a Whittaker differential equation is defined as [27] 

Two independent solutions of (42) are the Whittaker functions 

M«,^(x) = e-ix5$(a, c; x), (43) 
W^K,^(a;) = e-tarf *(a, c; x), (44) 

where 
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^ and $ are confluent hyper geometric functions. The solution of Eq.(38) that is regular at 
r — oo, has the form 



(47) 



where u — \J^i — \i which vanishes when Nf Nc- Eq.(39) has a solution of the type 
(43), which is regular at r = 0: 



^2(r) =e-'^''$(l-^,2;27firj , r < 1. 



(48) 



Since these two different solutions represent one solution of our original problem (21), we 
have to match them at some intermediate point r = ro = 0{1). A choice of ro corresponding 
to Appclquist et al. would be ro = 1/a — 1. But for the time being we keep an arbitrary 
matching point. As matching condition we take 



dr 



log 



0. 



r=ro 



With 



$'(a, c; x) = -$(a + 1, c + 1; x), 
c 

^''(a, c; x) = -a^'(a + 1, c + 1; x), 
the matching condition (49) reads 

/I As \ ^ (I - ^ + iz/, 2 + 2iz/; 2mro 

2m - — — + — 7-; 

V2 2 y * (1 - ^ + iz/, 1 + 2ii/; 2mro 

1 /I . \ ^~ (~ A.\ $(2-ir.3:2/m„ 
— ] — m + m — \ m — 



(49) 



(50) 
(51) 



ro V2 



2/ $(l-^.2;2mro) 



(52) 



Since we consider small m, the dependence on m on right-hand side of the previous equation 
vanishes. We define 



, . $(2-||,3;2mro) 

a — m—\m — — — ) ; (-, 

V 2y $(i_^,2;2mro) 



(53) 



with rfi = y ^ . Now 2mro -C 1 , so for the ^-function we can use the asymptotic form for 
small argument, i.e. 
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^(a, c; x) 



Til - c) 



■=-0 r(a-c + i; 

The matching condition then reads 



'1 A3^.^ 

2-1-+'", 



r (-1 - 2iu) 

.r(|-^--) 
r {-2w) 

which can be rewritten as 



r(c-i) 
+ r(a) 



r (1 + 2iz/) 



(2mro)^+^'^ + ' "^"^ X (2mro)- 



1 + 



r (2iz/) 



(2mro)'" + 



1 



3 _ 

2 2 

r {2w) 



r i-f + 



r (-2iz/) 



(2mro)' 



V2 2 



(2mro 



= dro 



With the use of 

ae'"^ + ae-''^ = 2\a\ cos(0 + 9), 9 = arg(a), 
the equation (56) can be reduced to the following one 



where 



E(z/) = arg 



' r(l + 2w) 



Eq.(58) then gives for m: 



m — exp 



2v 



1, 



- log(2ro) + - arctan — h -S(z/) 

V 1 — 2aro 

with n a positive integer. In the limit — > 0, Eq.(60) takes the form 

27rn 



m = exp 



where 



6 = - log(2ro) 



1 - 2dro 



27 + ^ 



1 A, 



with the di-gamma function, ■^(z) = r'(2;)/r(z), and where 7 is the Euler constant, 
the solution with the largest value of m < 1 {n—1) corresponds to the ground state, 
it has the lowest energy. 
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With the help of Mathematica we can calculate b explicitly for various values of the 
parameter a {tq = 1). With u — 



Ai = ^, A. = -, A3 = ^, m = y-, (63) 

this gives 

massless fermion loops: (7 = 0, b — 2.214, 

Stt 

massive fermion loops: <^ = b — 7.136. (64) 

Thus we have found that for a vacuum polarization with massive fermions, the infrared 
dynamical mass m obeys an exponential scaling law near the critical point N^. despite the 
presence of a Coulomb tail in the potential (27). This means that chiral symmetry breaking 
takes place essentially at intermediate distances 1 < r < 1/m. For the massless vacuum 
polarization, this was first pointed out by Appelquist ct al. [2]. The main difference between 
a vacuum polarization with massive (parameter cr 7^ 0) or with massless fermions (a = 0) 
lies in the scaling factor b, which is larger in the case of massive fermion loops in the vacuum 
polarization. We note here that our real space approach automatically takes into account the 
ultra violet tail of the integral equation in momentum space to which an attention was paid 
in Ref. [28]. Moreover, in comparison with their massless vacuum polarization calculations, 
our analytical solution gives the scaling factor b closer to exact computation of that from 
whole nonlinear equation (compare our b = 2.214 with Kondo and Maris' b = 3.94 [28] 
against exact b = 1.842). 



V. NUMERICAL CALCULATIONS IN MOMENTUM 

SPACE 



In this section we discuss the numerical calculation of the scaling law of the dynami- 
cal mass near criticality for a vacuum polarization with respectively massive and massless 
fermion loops. The nonlinear DSE can be written in the following way, with a = 1 



dk — . L.. „s K{p\ k\ m). 



(65) 



where the kernel K is given by 



K{p\k\m) = / d^ sin^ 



1 



l + n(g2) 2 Jo l + U{q^y) 



dy 



(66) 



with q"^ = k"^ + — 2pk cos 9, and we recall the expression for the one- loop vacuum polar- 
ization. 



2m + 



f — 4m^ 



arctan ( ) 
q \2mJ 



(67) 
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In the case of massless fcrmions in the vacuum polarization, the y-integration in (66) can 
be performed exphcitly. It reduces to 



dO sin6'— 

r 



l+i_J_log(l + ,) 



(68) 



which docs not depend on the dynamical mass, so that it only has to be calculated once. It 
is convenient to perform the calculations on a logp^ scale in order to see more details of the 
structure. 



M{x) 



^ logA2 

log//2 



y + M^{y) 



K{x,y,m), 



(69) 



where x = p'^, and y = k"^ = exp(t). We have introduced infrared and ultraviolet cutoffs, // 
and A, respectively, which should satisfy <C (0), and A 3> a = 1. 

The basic idea of our numerical calculations is that we solve Eq.(65) iteratively as an 
integral equation. First we calculate the kernel K for all grid points x and y with some initial 
value for m. Then we iterate the integral equation until it has reached some convergence 
criterion. Then we recalculate the kernel with the new m = M(0), and iterate again. This 
procedure we repeat to the point were the complete system has converged. To obtain a 
satisfactory numerical procedure, the integrals in K are performed using Gauss-Legendre 
quadrature, which works well since our integrands in (69) are sufficiently smooth. For M 
at the grid points x we use cubic spline interpolation on a logx-scale. The ^-integration is 
performed using Gauss-Legendre quadrature between the grid points x, 



M{xi) 



/ * 

log Xj 



yyyMjy) 

y + M^{y) 



K{xi,y,m) 



(70) 



where xi — /i^, and Xn — A^. With this procedure we find that for calculation near the 
critical point Nc, we need at least 128 grid points for x. We have calculated the mass 
function M{x) for a series of values of Nf near A^^- For all cases it appears that the infrared 
m obeys the exponential scaling law. This has been verified using a least squares fit for the 
numerical data. We have fitted the infrared mass m, as a function of Nf, with the following 
function 



logm(A^/) = 



27ra 



ic/Nf - 1)' 



+ b, 



(71) 



where a, b, c, and d are the parameters to be fitted. For the nonlinear massless case the 
least squares fitting of the data gives 



a = 0.9907, 6= 1.720, c = 4.32312, d = 0.5018, 
and with massive fermion loops 

a = 1.055, b = 3.494, c = 4.32315, d = 0.4893, 



(72) 



(73) 
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which is shown in Fig. 2. This fitting proves the exponential scahng law (61). To obtain 
values of b and compare them with the analytic calculation we take a — 1, c — Nf. — 128/37r^, 
d — 1/2. The numerical results are: 

massless fermion loops: b — 1.842, 

massive fermion loops: b — 2.955. (74) 

We also investigated the difference between a linearized approximation (18) and the nonlin- 
ear DSE (15). Qualitively the linearized approximation is equivalent to the more realistic 
nonlinear equation. Quantitively the scaling factor 6 is a few percent smaller for the lin- 
earized model. 

In Fig. 3 we have plotted the dynamical fermion mass function M{p^), computed for 
various values of the fermion number Nf. Apart of Nc it has behavior reminiscent to that in 
the quenched approximation [9,10]: the mass function is constant at momenta p <^ a and 
decreases as o? j]? at large momenta. The mass function changes its behavior at momentum 
p ~ a, a natural scale of the model. Close to critical number of fermions, N j ^ A^cj we 
observe three different regions in behavior of M(p^): M{p^) ^ m at p < m, M{p'^) ~ 
vfV'l'^ j y/p oXm <p <(x and M(p'^) ~ m^/p^ at p > a. The scale m is exponentially small in 
comparison with the intrinsic scale a, see Eq.(71). Fig. 3 shows convincingly how a hierarchy 
of scales can occur in a theory under specific conditions (in our case when Nf is close to Nc) . 
This might be important for further understanding a hierarchy of scales in unified theories. 

VI. CONCLUSION AND DISCUSSION 

In the present paper we investigated xSB in (2-|-l)-dimensional QED with Nf number of 
four-component Dirac fermions. We have pointed out one possible scheme for truncating an 
infinite chain of SDE, which is consistent with WT identity and which seems to have been 
overlooked so far. Further simplification of that coupled system of SDE led to studying the 
SD equation for a fermion mass function where the full vacmim polarization with dynamically 
massive fermions has been taken into account. The latter proves to be important for keeping 
the confinement property of the model at large distances. The study of the linearized 
equation for a fermion mass function has been performed in real space, and that allowed us 
to avoid making drastic approximations in order to get analytical results. Our analytical 
solution fits nicely numerical calculations of the full nonlinear integral equation for the mass 
function. 

Further investigation of aforementioned truncated system of SDE is now in progress. 
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FIGURES 

FIG. 1. Truncated system of Schwingcr-Dyson equations. 

FIG. 2. The numerical results for M (0) as function of Nf for massive and massless fermion 
loops in vacuum polarization. 

FIG. 3. Dynamical mass function M{p'^) for Nf = 1.0, 2.0, 3.0, 4.0. 
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